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ABSTRACT 

Centaurus (Cen) A represents one of the best candidates for an isolated, compact, highly polarized source that 
is bright at typical cosmic microwave background (CMB) experiment frequencies. We present measurements 
of the 4° x 2° region centered on Cen A with QUaD, a CMB polarimeter whose absolute polarization angle is 
known to 0°5. Simulations are performed to assess the effect of misestimation of the instrumental parameters 
on the final measurement, and systematic errors due to the field's background structure and temporal variability 
from Cen As nuclear region are determined. The total (g, U) of the inner lobe region is (1.00 ± 0.07 (stat.) ± 
0.04 (sys.), -1.72 ± 0.06 ± 0.05) Jy at 100 GHz and (0.80 ± 0.06 ± 0.06,-1.40 ± 0.07 ± 0.08) Jy at 150 GHz, 
leading to polarization angles and total errors of -30?0 ± 1?1 and -29? 1 ± 1°7. These measurements will allow 
the use of Cen A as a polarized calibration source for future millimeter experiments. 

Subject headings: cosmic background radiation - galaxies: individual (Centaurus A) - instrumentation: po- 
larimeter s - radio continuum: general 



1. INTRODUCTION 

Cosmic microwave background (CMB) polarimeters re- 
quire accurate calibration of both their temperature and polar- 
ization properties. Unfortunately, astronomical sources suit- 
able for polarization calibration in the wavelength range at 
which CMB experiments typically operate are rare. With a 
number of exquisitely sensitive instruments in the field or near 
to deployment, the need for astrophysical polarization calibra- 
tion sources useful for instruments with beam sizes of several 
to tens of arcmin is acute. 

An ideal polarized calibration source would be compact, 
bright in both temperature and polarization, static over time, 
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and isolated from galactic diffuse emission. Unfortunately, 
the typical frequencies at which CMB experiments are de- 
signed to measure are minima in the emission spectra for 
other astrophysical sources of radiation (synchrotron, free- 
free, and thermal dust emission), so such sources have not 
been identified in large numbers or studied in great detail at 
the frequencies of interest. A noteworthy candidate calibra- 
tion source which has been observed with CMB polarimeters 
in the past is the superno va remnant Taurus A (Leitcl TeTaLl 
2002, Barkats et a 1. 2005), though as an object in the northern 
hemisphere it is not available to telescopes at all sites. 

One of the best candidates for a CMB polarization cali- 
bration source in the southern hemisph ere is the ra dio bright 
galaxy Centaurus A (hereafter Cen A; llsraell fl998l provides 
an excellent review of its multi- wavelength properties). The 
optical counterpart of Cen A is the massive elliptical galaxy 
NGC 5128 at a distance of 3.4Mpc. At radio frequencies Cen 
A presents rich structure over many decades in angular size 
driven by its nuclear source. Collimated radio jets are emitted 
from the compact nucleus and become sub- sonic a few parsec 
(~ 0.1") from the central source. At 5 kpc (~ 4') from the 
nucleus, the jets expand into plumes which spread up to 250 
kpc (~ 3°) into the inter cluster medium (ICM). The interface 
between the bright sub-sonic jets and the smooth, low sur- 
face brightness plumes is abrupt, allowing high signal to noise 
measurements of the jets themselves. The ~ 10' jets, known 
as Cen As 'inner lobes' in the literature, have a spectral index 
of -0.7 from > 500 MHz to ~ 100GHz and are both polarized 
and bright at v > 100 GHz. Moreover, Cen As inner lobes are 
unique as their properties do not vary on human time scales, 
which makes them perfect candidates for polarization calibra- 
tion measurements. Cen A also lies well above the galactic 
plane, is well matched to the size of typical CMB polarime- 
ter beam sizes and is well studied at many other wavelengths, 
particularly in the radio. In this paper we present the results 
of observations of Cen A with QUaD, a CMB polarimeter ca- 
pable of simultaneously measuring Stokes' /, Q and U param- 
eters with bands at 100 and 150 GHz and angular resolution 
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5'.0 and 3 '.5 respectively. 

2. INSTRUMENT SUMMARY AND OBSERVATIONS 

A full descript i on of the QUaD instrument can be found 
in lHinderks et al.1 (120091) (hereafter H09). QUaD was a 2.6m 
Cassegrain telescope using the pre-existing D ASI mount at 
the US south pole station (iLeitch et alJl2002h . QUaD oper- 
ated from 2005 February to 2007 November, taking astro- 
nomical data during the austral autumn through winter into 
spring; maintenance and calibrations were performed during 
the austral summer. The receiver comprised 3 1 pairs of polar- 
ization sensitive bolometers (PSBs), 12 pairs at 100 GHz and 
19 pairs at 150 GHz. At either frequency these were split into 
two polarization groups which measure either Q or U\ each 
bolometer pair is also sensitive to Stokes' /. 

The data presented here were obtained on 2006 June 7 and 
2007 October 2, 3 and 4. The observations were performed 
in a similar m anner to a standard QUaD observation block 
as described in iPryke et al.l (l2009h (hereafter P09), although 
the lead-trail scan strategy used in other QUaD observations 
is not employed. Each scan group begins with an elevation 
nod during which the telescope is moved up and down by 
1° in elevation to inject a signal of common amplitude into 
each bolometer. The telescope is then scanned back and forth 
at constant elevation over a throw of 3?2 in sidereal tracking 
corrected azimuth; four such scans constitute a single scan 
block. For these observations, the scan rate is 0°066 per sec- 
ond in azimuth, or 0?05 per second on the sky at Cen As ele- 
vation. After each constant elevation scan block, the elevation 
is stepped by 1.2 arcmin and the next block is begun. In all, 
101 scan blocks were performed, equating to a rastered data 
set mapping ~ 4° x 2° on the sky to a uniform noise level. 
These 101 scans required 15.6 hr; adding the standard QUaD 
calibration set before and after the scan block (see H09 for a 
description of such a set) leads to a total observation time of 
17.6 hr per day. 

The DASI mount supports rotation of the entire telescope 
about its boresight; on a given day this rotation was fixed for 
the Cen A observations. However, the rotation was varied 
over different days, which has the benefit of changing the po- 
larization characteristics of the instrument in a known manner. 
The angles chosen were {0°,-30 o ,+30 o ,0°} for 2006 June 7 
and 2007 October 2, 3, and 4, respectively. 

3. FROM TIME STREAMS TO T MAPS 

Unlike standard CMB observations, the QUaD polarization 
maps of Cen A require systematic effect corrections derived 
from the temperature maps (Section 0]), so it is necessary to 
utilize a two step process where the total intensity maps and 
then the polarization maps are constructed. The Cen A tem- 
perature maps themselves can be constructed directly from the 
data time series without any such corrections. 

In order to produce maps, the relationship between the mea- 
sured time series dit) and the sky temperatures T,Q,U must 
be known. Neglecting noise contributions and the bolometer 
transfer function, it can be written as 



d(t) = g(t) / drB(r-r b )x 



1-e. 



T(r) + Y^r e i 2(r) cos 2^ + U(r) sin 2^} 



(1) 



where r is the direction of observation and ip is the detec- 
tor's polarization angle (the time dependence of these two 



quantities is suppressed for clarity). The time dependent gain 
g, cross polar leakage e and polarized beam shape B are all 
quantities which must be determined via calibration measure- 
ments, as must be the zero point of ip. 

Initial low-level processing of the time series is performed 
using the same algorithm as presented in P09. The bolome- 
ter time constants and electronic filters are first deconvolved 
to reverse their effect on the time series. These time series 
are then despiked by removing those individual scans which 
exhibit transient impulse events from the analysis. The eleva- 
tion nods are used to determine g(t) for each scan block by us- 
ing the changing atmospheric loading to compute the relative 
responsivity of each bolometer in a frequency group. Each 
detector in a frequency group is scaled to the group's mean, 
yielding relatively calibrated time series for the entire array. 

Three more pieces of information are required to create 
maps from the rectified and responsivity calibrated bolome- 
ter time ordered data (TOD): the telescope pointing, the PSB 
angles, and the PSB efficiencies. Knowing the pointing of the 
telescope requires knowledge of both the sky offset of each of 
the detectors and the absolute pointing of the boresight, the 
derivation of these parameters is described in detail in H09. 
The detector offset angles used in the map making process are 
derived from monthly measurements and have an estimated 
uncertainty of 0M5; since the signal-to-noise ratio on Cen A 
is large the absolute per day pointing solution can be directly 
determined from / maps. 

The PSB polarization angles and efficiencies were deter- 
mined using a chopped thermal source placed behind a polar- 
izing grid on a mast near the telescope; H09 reviews the mea- 
surement and preliminary analysis of these data. The cross 
polar leakage is defined as the ratio of response to anti-aligned 
and co-aligned incident polarized light for a given PSB. The 
measured values of e have a mean of 0.08 with an rms scatter 
of0.015. 

Knowledge of the error on the absolute polarization angle 
has improved over that presented in H09 thanks to both an 
improvement in the analysis of the calibration data and a new 
method of measuring this angle. The new method of angle 
measurement relies on the fact that the polarized response of 
the detectors is mechanically constrained to lie along lines of 
symmetry of the detector rows, any misalignment is estimated 
to be random and < 0°.5 per detector. Measurement of the 
polarization angle is performed by scanning rows of detectors 
across bright sources at telescope rotation angles where the 
polarized response of the detectors matches the orientation of 
the rows. As there are several such symmetries, a number 
of independent measurements can be made. The final error 
ascribed to this angle measurement method is 0?15. 

The improved analysis of the calibration source data pre- 
sented in H09 is based on the observation that the orientation 
of the polarized source was not well constrained in the az- 
imuthal plane; azimuthal misalignments at the source box can 
result in spurious polarization measured at the telescope. We 
therefore discard all measurements performed with the cali- 
bration source grid in the horizontal orientation and use only 
the vertical grid measurements; these are immune to such 
alignment errors and have a well measured orientation with 
respect to gravity. The updated analysis which uses only the 
grid- vertical measurements has a scatter of 0?2 over the mea- 
surements. The absolute polarization angle of the telescope 
derived from these measurements is in good agreement with 
the angle derived from the row scans. We therefore use this 
revised absolute telescope angle and conservatively ascribe an 
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error of ±0?5 to its measurement. 

To construct maps further processing of the TOD is re- 
quired. The PSB pair time streams are first summed to cre- 
ate total intensity TOD corresponding to the Stokes' param- 
eter / and differenced to make instrument frame Stokes' Q 
or U ; these must be rotated to an absolute reference frame 
for polarization using the known polarization angle of the 
telescope (as discussed in Section [4]). Drifts in the atmo- 
spheric emission and instrument create 1/f noise in the TOD 
much larger than photon noise on long time scales; to re- 
move these correlated signals a fifth-order polynomial is fit 
to each 48 s half- scan and subtracted from the sum and dif- 
ference TOD. This is a higher order than that chosen for the 
CMB data since we are not attempting to measure large spa- 
tial modes in the Cen A maps. As the signal-to-noise ra- 
tio on the source is large in a single scan, it is necessary to 
mask source flux to prevent biasing the polynomial fits. The 
mask used here comprises a circular region with r = 0°.2 cen- 
tered on Cen A; this size was chosen to preserve structure 
in the inner lobe region while maximally removing drifts in 
the TOD. The same mask is used for the sum and difference 
TOD. When filtering the sum data to produce / maps, the 
fit mask is also augmented by two circular regions centered 
at R.A.= {13*450, 13*353}, decl.= {-42?80,-43?70}; these 
mask Cen As moderately bright northern spur region and the 
elliptical galaxy NGC 5090 from the polynomial filter. 

Maps are made by binning the TOD into a grid of pixels 
weighting by the inverse variance of each half- scan computed 
after the polynomial filter has been applied. The map pix- 
elization is in R.A. and decl. using square pixels 0°02 on a 
side. The per day signal-to-noise ratio in the / maps is large, 
so by comparing individual day maps it is possible to observe 
the telescope pointing solution vary by up to 30 arcsec day to 
day. To correct for this variation, the per day / maps are made 
and that day's astrometric offsets are found by minimizing 
the sum of the squared pixel differences in the region of Cen 
A between that day's map and the map derived from the first 
day's observation with its astrometric calibration manually set 
to Cen As known position. Once determined, these offsets are 
then always applied to both sum and difference TOD during 
the map making process. 

To calibrate to CMB mK we utilize the CMB power spec- 
trum c ross calibration procedure presented in iBrown et"aTI 
( 2009|); the resulting temperature maps are shown in Figure 
[U To calibrate between / maps in Jy sr -1 and T maps in ther- 
modynamic temperature the relation 

dI = (dB/dT) 2J3K dT (2) 

is used. These values correspond to {29.2,58.5} Jy map 
pixel" 1 K- 1 at {100, 150} GHz. 

4. POLARIZED INSTRUMENTAL EFFECT CORRECTION 

There are a wide range of systematic effects which can mix 
/ to polarization and Q to/from U in polarimeters. The QUaD 
simulator, discussed in detail in P09, is an extremely detailed 
model of the instrument and various sources of noise whose 
primary purpose is to model time s treams for use in M ASTER- 
style CMB power spectra analysis dHivon et al.ll2002|) . As the 
simulator is built to model the instrument accurately it can 
be used to measure and correct for instrumental systematic 
effects which are known to exist in QUaD data, whatever the 
astronomical source under study. In order to produce clean 
polarization maps, it is necessary to simulate and correct for 
these polarization mixing effects. 



4.1. Polarization leakage correction 

Although the detector offset angles show no evidence for 
variation over time or between bolometer pairs sharing a 
feed horn, they do exhibit repeatable offsets between the two 
halves of each detector pair with an rms magnitude of 0M 
over the array. When the sum and difference maps are con- 
structed, the mean detector offset of a pair is used; unfortu- 
nately, the difference between this mean and the actual point- 
ing of a detector can mix / — >> Q,U. This polarization leakage 
can be quantified using simulations of the observations. 

To measure polarization leakage, an input map is generated 
from the coadded / maps at each frequency shown in Figure 
[U Hereafter, we denote the "polynomial filter mask region" to 
be that area in the map corresponding to the regions masked in 
the time stream polynomial fit discussed in Section i.e. the 
r < 0°.2 area centered on the source. Regions outside of the 
polynomial filter mask region are set to zero, such that only 
the central area contains non-zero flux. Each detector's point- 
ing time stream is constructed using the individual measured 
offsets, which vary between detectors in a pair. These time 
series are used to sample a new, simulated time stream from 
the measured / maps. As no noise is added to the TOD and the 
polarization maps are set to zero, this procedure generates a 
simulation of the pure / component of the Cen A observation 
for each detector. The PSB pairs are then differenced in the 
standard way; if the pointing between pairs was perfect, this 
procedure would cancel the structure appearing in the con- 
structed / maps perfectly. However, as the pair pointing is 
not identical, this procedure produces residuals reflecting the 
leakage of / to polarization in the actual observation. These 
leakage time streams are subtracted from their corresponding 
pairs in the data TOD, which corrects for the leakage effects 
at the time sample level. 

These polarization leakage corrected time streams are 
mapped onto the same R.A. and decl. grid as the / data. For 
difference time stream data, the product of the data and the 
sine and cosine of the detector angles as projected on the sky 
must be accumulated for each pixel to weight the different 
detector pairs according to their polarization sensitivities; this 
procedure rotates instrument frame Q and U to absolute Q 
and U. The set of 2 x 2 matrices comprising products of these 
sines, cosines and difference time streams must be inverted for 
each pixel; these are then accumulated into the grid to pro- 
duce absolutely referenced Q and U maps. In this work we 
follow the IAU polarization convention with positive Q run- 
ning north-south and positiv e U running northeast-southwest 
(Ha maker & Bregmanl ll996). Figure [2 shows the CMB tem- 
perature calibrated Q and U maps at 100 and 150 GHz. These, 
combined with the / maps, are the fundamental output of the 
analysis pipeline. 

These maps can be used to characterize the polarization of 
Cen A per QUaD beam or per map pixel as desired (see Sec- 
tion [5]). 

4.2. Polynomial filter correction 

A second correction which needs to be applied to these data 
is due to the time series polynomial filter. The filter mask dis- 
cussed in Section [3] is constructed as a compromise between 
maximizing the mask size to exclude as much source flux as 
possible from the polynomial filter fit and the need to mini- 
mize the mask so that individual scans are well constrained 
when the polynomial function is determined from the data. 
Although the mask effectively excludes the bright source flux 
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Fig. 1 . — Temperature maps of Centaurus A. Temperature scale is /xKcmb- Top: 10 GHz (left) a nd 150 GHz (righ t) maps of the Cen A region at low contrast. 
The R.A. and decl. range has been chosen to facilitate comparison with Figure 3 of Israel (1998) and Figure 1 of Israel et al. ( 2008); none of the large scale 
structure visible in those maps at lower frequencies is present in the QUaD data. The structure of the inner lobe region comprising two synchrotron emitting 
convexities in the ICM is unresolved at 100 GHz but resolved at 150 GHz. Bottom: 100 GHz (left) and 150 GHz (right) maps of Cen A at high contrast. 
Low surface brightness emission is visible to the north east of the inner lobe region, as is NCG 5090, an unrelated elliptical galaxy in the Cen A field near 
aj2000 = 200°266, 5 moo = -43°741. 



during the polynomial filter fit, there always remains some 
small fraction of the source flux outside the mask for the scan 
length used in these measurements. This flux has the effect of 
biasing the filter determined from the data; this can be quan- 
tified using the QUaD simulator. 

In order to produce realistic simulations of these observa- 
tions, it is necessary to create an accurate source map at each 
polarization and frequency. Such a map needs to have high 
resolution so that the effect of the QUaD beams can be accu- 
rately simulated, and must have large enough signal-to-noise 
ratio on all components of the source that the input map is 
injecting negligible error into the simulation. In the case of 
the CMB highly accurate simulations of the background sky 
based on known physical mechanisms at the epoch of recom- 
bination exist; unfortunately, the same is not true of Cen A. 
Such a map of the Cen A field at QUaD's frequencies can 
be produced either by scaling a map from a different fre- 
quency by the relevant spectral index of the source, or by 
using the QUaD maps themselves. The former method is dif- 
ficult: maps of the Cen A region at radio frequencies tend to 
either have much higher resolution but are not sensitive to ex- 
tended emission, or much lower resolution which would not 



allow us to resolve the source. We therefore have utilized the 
latter approach. 

In order to increase the angular resolution of the Cen A 
maps, an algorithm which deconvolves the effective beam 
from the observed maps is necessary. We have chosen the 
Richardson-Lucy deconvolution (RLD) algorithm here as it 
can efficiently deconvolve a k nown kernel from noisy data 
(iRichardsonll 19721: iLucvll 19741) . The RLD algorithm employs 
B ayes' theorem to iteratively reconstruct the maximum likeli- 
hood background map given the beam kernel and the observed 
map with which it has been convolved. For the QUaD mea- 
surements of Cen A, the input maps are the observed 4 day 
coadded maps shown in Figures Q] and [2] 

The effective QUaD beam for these observations has been 
calculated using noiseless simulations whose input is a delta 
function. These computed beams are well matched to a Gaus- 
sian to a level below -15 dB, so for computational ease we 
choose a RLD beam kernel to be a symmetric Gaussian with 
FWHM matching the QUaD beams at either frequency (H09). 
Applying the RLD algorithm to the 4 day coadded Cen A 
images yields the maps shown in Figure The deconvolu- 
tion is performed for each polarization state at both frequen- 
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Fig. 2. — Polarization leakage corrected maps of Centaurus A. Temperature scale is /xKcmb- Top: 100 GHz (left) and 150 GHz (right) maps of the Cen A 
region for Stokes' Q. Bottom: 100 GHz (left) and 150 GHz (right) maps of the Cen A region for Stokes' U. The polarization of the northern inner lobe is about 
20 % of / at either frequency. The polarization fraction in the southern lobe is much smaller. 

cies; in the case of Q and U the polarization leakage corrected 
maps are used in the deconvolution so that polarization leak- 
age need not be simulated. 

The polynomial filter correction is determined by perform- 
ing simulations similar to those discussed in the previous sec- 
tion. The RLD input maps are convolved with the ideal instru- 
ment response function at each frequency and TOD are con- 
structed for each detector using the known telescope pointing. 
These TOD are then summed and differenced pair- wise, and 
a fifth order polynomial is fit and subtracted from each pair's 
time streams. These filtered TOD are then binned in the usual 
way to produce fully corrected polarization maps; these con- 
stitute QUaD's final processed maps of Cen A. 



A quantity which encapsulates a great deal about these po- 
larization maps is the total Q and U in an aperture. Here we 
employ a circular aperture with r = 0°.2 centered on Cen A's 
position; this matches the source mask used during the time 
stream polynomial filtering. The total polarized temperatures 
for the final 4 day co-added maps are given by 



V = {^Q(a,5)^U(a,5)}. 



(3) 



5. RESULTS 

Figure @] shows the 4 day co-added maps after the system- 
atic corrections have been applied. Using the RLD algorithm, 
these can be deconvolved with the nominal QUaD beam to 
remove the effect of the instrument transfer function to yield 
the maps similar to those shown in Figure [3} these maps al- 
low simulations of this source in /, Q and U for any instrument 
with beam size similar to QUaD's or larger 1 . 

1 Both systematic corrected as-measured and RLD maps, including esti- 
mates of the noise in each map pixel, will be publicly available on the web at 



The values of the aperture sums V measured from the QUaD 
maps are listed in Table Q] The error estimate from the map, 
denoted cr map CP), is computed by calculating the rms of pix- 
els in the same circular aperture offset from the source by 
AR.A. = 0?5 and dividing by the square root of the number 
of map pixels in the sum region. The error estimates cr map (V) 
are computed separately for Q and U at both frequencies. 

An independent measurement of the error in V can be pro- 
vided by the QUaD simulator, which includes an extremely 
accurate model of the noise in QUaD data. As in previous 
simulations, the RLD input map is convolved with the individ- 
ual detector's beams and sampled with the telescope pointing. 
Based on each detector's statistical properties over a noise 
measurement block, noise is generated and injected into the 

http : / / f ind. uchicago . edu/quad/quad_CenA/ . 
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Fig. 3. — RLD maps of Centaurus A. The deconvolution kernel is a sym- 
metric Gaussian with an FWHM matching the measured QUaD beam. Rows 
correspond to T (top), Q (middle), and U (bottom) while columns are for 100 
GHz (left) and 150 GHz (right). Temperature scale is /xKcmb; the scaling on 
the U maps has been reversed to allow comparison with the Q maps. Note 
that the deconvolved maps resolve the source into two distinct lobes with lit- 
tle or no polarized emission arising from the nuclear region between them, 
particularly at 150 GHz. 



TABLE 1 

Polarized Flux of Cen A from Corrected QUaD Q and U Maps. 





100 GHz 


150 GHz 




1.00 Jy 


0.80 Jy 


C r map(J]2) 


0.06 Jy 


0.05 Jy 




0.06 Jy 


0.06 Jy 




-1.72 Jy 


-1.40 Jy 


CTmap(E C/ ) 


0.08 Jy 


0.05 Jy 


&sim(J2 U) 


0.06 Jy 


0.07 Jy 





-30°0 


-29° 1 


CT(0) 


0°9 


1°2 



scan TOD. P09 details the noise model construction in detail; 
it is useful to note here that the noise estimate for each de- 
tector is constructed from statistical blocks of that detector's 
actual time series 5 scans long. These time streams are then 
fifth order polynomial filtered and binned into I,Q,U maps as 
with the observed data. This simulation process is repeated 
N times to obtain N random realizations of the measurement. 
The scatter of V in these N realizations yields the variation in 
the possible measurement outcomes due to the random noise 
in the data. 

Figure [5] shows the results of 256 such realizations for each 
day of observation. The standard deviations of these distri- 
butions are consistent with the aperture rms of the individual 
day maps, showing that the noise model is in good agreement 
with the variance in the data. Also plotted in this figure are the 
corrected per day V Cen A measurements; the hypothesis that 
these measurements are drawn from the same distribution as 



the simulated sample can be checked using Student's £-test. 
In the sample of 16 different permutations of frequency, Q, 
U and day, the null hypothesis is never rejected. This result 
is consistent with the null hypothesis for the overall set, and 
shows both that we have not detected a bias in the individual 
days' data and that the variance in the simulations is a realistic 
model for the variance in the real data. The random error es- 
timate on V derived from the scatter of the simulation results 
is listed in Table Q] as cr S i m (V). 

The overall polarization angle of the source can be com- 
puted using 

= ^arctan(^£//^e). (4) 

The values obtained in the QUaD measurement are listed in 
Table [TJ including the random error a(0). The released data 
maps and equation [4] can be applied per map pixel to deter- 
mine the polarization angle variation over the source. 

6. SYSTEMATIC ISSUES 

There are a number of systematic effects which can add 
error to the measurement of V and 6. These broadly group 
into either astrophysical effects or instrumental effects; these 
are discussed below. 

6. 1 . Background Structure 

The background structure in the Cen A field can potentially 
cause significant systematic error in this measurement. A 
number of measurements have mapped out Cen As 8° x 4° ra- 
dio emitting outer lobes at various freque ncies, most recentl y 
at high rad i o freq uencies with WMAP (I Wright et al.l l2009|) . 
Ilsrael et al.l (120081) show detailed WMAP maps of Cen A at 23, 
33, 41, 61 and 94 GHz; the outer lobe emission is clearly visi- 
ble at the lower frequencies but is not significantly detected at 
94 GHz. This result is evidence that at the QUaD frequencies 
Cen As outer lobe synchrotron emission is much below the 
brightness of the inner lobe region, but this result should be 
checked. 

To perform this assessment, maps are first made in the stan- 
dard way, except a first order polynomial filter is used to re- 
move the instrumental drifts from the time series of each de- 
tector; unlike the standard fifth-order polynomial used in the 
normal map making process, this filter preserves structures 
larger than 1 ° on the sky. The central r < 0°.2 region is again 
masked during the polynomial fitting procedure to avoid bias. 
These data are binned into maps, and then the central poly- 
nomial mask region is replaced with the mean of the pixels in 
the annulus 0?2 < r < 0?3 just beyond the mask region. These 
new maps are then convolved with a Gaussian kernel with an 
FWHM of 1 ° ; this smooths out the small features whilst re- 
taining the largest structures in the map. Figure [6] shows the 
resulting maps for T, Q and U. Of these, only the 100 GHz 
T map shows evidence for background structure on scales 
~ 1°. The total flux in the region where T > 20/iK in this 
background-filtered 100 GHz map is 10% of the total flux in 
the inner lobe region at 100 GHz. A number of different ker- 
nel sizes from 0°.2 to 2° have been applied in this procedure 
and the total flux in the background region does not change 
appreciably under different kernel widths. Table [2] lists the 
estimated uncertainty in V and caused by the background 
structure; for each polarization and frequency this quantity is 
calculated by summing the background signal in the region 
where Tiooghz > 20/iK and multiplying the total by the ratio 
of the area of the r < 0?2 region to the area of the summed 
region. 
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FIG. 4. — Final T maps with fully processed polarization vectors over plotted. The left- and right-hand panels show the 100 and 150 GHz maps, respectively. 
Polarization vectors are plotted for those map pixels with T > 1000/iK. The color scale is in /iK and the length of a 4mK pure +U vector is indicated. 
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Fig. 5. — Results of 256 simulations of the per day Cen A observations for 
Q and U at 100 and 150 GHz. Days run row-wise and polarization and fre- 
quencies run column- wise. The simulations use knowledge of the individual 
detectors' beams, the full QUaD noise model, and share the same analysis 
pipeline as the actual data. The individual realizations of V at both frequen- 
cies are shown in the histograms; also shown are the input V (dotted lines) 
and the V measured from the data (solid lines) for each day, frequency, and 
polarization. This result shows that there is no bias in the simulations, that 
the measured data are compatible with random draws from such realizations, 
and that the noise in V can be estimated by the scatter in the simulation real- 
izations. 

6.2. Source Variability 

Another potential source of systematic error is temporal 
flux variability in Cen A itself. Although the inner lobe re- 
gion is too large for its emission to change on human time 
scales, the nucleus of the source is known to be variable in /. 
Further, no measurements of this source's variability in polar- 
ization around 100 GHz exist. 

To measure the /, Q and U variability in Cen A, the nomi- 
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Fig. 6. — Background-filtered maps of the Cen A region. The left-hand 
column shows 100 GHz maps for T (top), Q (middle), and U (bottom); the 
right-hand column is the same for the 150 GHz maps. Color scale is in /iK 
and the cross in each map shows the nominal position of Cen A. The 100 GHz 
T map is the only one exhibiting evidence of background structure; there is 
no evidence for background contamination in the polarization maps. 

nal QUaD beam shape centered on the Cen A nucleus is fit to 
the final processed beam convolved maps at each frequency. 
As the size of Cen As nucleus is much smaller than a QUaD 
beam, this fit yields the flux in the nuclear region, and sub- 
tracting it from the map yields the total flux in the inner lobe 
region for each frequency and polarization. The error in the 
resulting nuclear flux is calculated by measuring the rms de- 
viations of the pixels in a number of circular apertures which 
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have the same angular extent as a QUaD beam. The set of 
aperture rms statistics is then averaged to yield the noise in 
the measurement of the flux in a single aperture. Figure [7] 
shows the / measured by QUaD in Cen A's nuclear region 
for each day of observation at 100 and 150 GHz. At both 
QUaD's bands, the / flux from Cen A's nucleus is ~ 35% of 
the total flux from the inner lobe and nuclear regions; this im- 
plies that the temporal variability expected in the QUaD maps 
is 5, 10% in I. Also plotted are the nuclear fluxes measured 
by llsrael et al.l (l2008[) . These measurements were performed 
with a heterodyne system on SEST at several frequencies be- 
tween 80 and 300 GHz: in order to create data comparable 
to the QUaD set, we compute the QUaD bandpass-weighted 
average of the SEST points. 
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Fig. 7. — Variability in Centaurus A's nucleus at 100 and 15 GHz for 
Stokes' /. The open circles show data presented in Israel et al. ( 2008) for 
the period 1988-2004 obtained using a heterodyne system on the SEST tele- 
scope, while the dots show the QUaD data presented here. These data are 
compatible with the well-known ~ 30 % peak to peak radio variation of this 
source at these wavelengths. 



The llsrael etail (l2008h points clearly exhibit the w 30% 
peak to peak variation in / well known from measurements 
at lower radio frequencies. The QUaD measurements also ex- 
hibit these fluctuations, although essentially only at two inde- 
pendent times as the 2007 measurements were performed on 

3 consecutive days. These data do place a limit on the short 
term (day time scale) / variability of the Cen A nucleus of 1 % 
at 100 GHz and 6% at 150 GHz. The / flux of only the in- 
ner lobes in the standard r < 0°.2 aperture is measured to be 
14.2 ± 0.3 Jy at 100 GHz and 16.6 ± 1.1 Jy at 150 GHz. 

Cen A's nuclear flux in Q and U is computed in the same 
way as for /; Figure [8] shows the resulting Q and U for the 

4 days of observation. Comparing to the model that the flux 
of these 16 measurements is constant, the obtained reduced 
^ 2 = 0.64; based on the QUaD data, there is no evidence of 
polarized temporal variability in Cen A. An upper limit on the 
possible variation can be obtained by computing the mean ab- 
solute deviation of these individual day data from their mode; 
the limits are listed in Table [2 These limits are significantly 
smaller than the estimated 30% variation expected from the 
/ variability. Inspection of the RLD maps in Figure [3] shows 
that very little of the polarized emission is originating from 
the nuclear region itself. This is evidence that this source's 



polarization stability arises because the polarized millimeter- 
wave flux is predominantly being emitted by the inner lobe 
region rather than the nucleus. 
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Fig. 8.— Time variability of Cen A in total Q and U at 100 and 150 GHz 
(left and right, respectively). Symbol "+" shows Q and symbol "x" shows 
U for each of the four days of measurement; flux is given in Jy. The mean 
of each frequency and polarization group has been subtracted to isolate the 
variable component of the flux. As these points are consistent with no time 
variation, half of their peak to peak scatter is used as the error associated with 
variability in the measurement of Cen As polarization properties. 



6.3. Instrumental Parameters 

There are a number of instrumental parameters used in the 
QUaD analysis pipeline which can cause systematic errors; a 
review of these is given in P09. Following that work, for these 
Cen A data the possible systematic errors caused by misesti- 
mates of these parameters are investigated using the QUaD 
simulator. 

The absolute calibration unce rtainty from detect or output 
to Tcmb is estimated to be 3.4% (Br own et al.ll2009|) . As this 
uncertainty is completely correlated between /, Q and U it 
does not affect polarization angle or fraction, only absolute 
polarization temperature. For this reason it is omitted from 
Table O though formally it should be included in a complete 
accounting of the uncertai nty i n Cen A's Q and U. 

As discussed in Section 14.11 the detector pair beam offsets 
cause leakage of / to polarization. Although we have ac- 
counted for this effect, the measurement of the offsets is still 
subject to error. To measure the effect of a random misesti- 
mation of the detector offsets, we have performed simulations 
where the measured detector offsets are randomized by the 
rms centroid uncertainty of 0M5. Though this leads to no 
change in the mean value of V, it increases the noise in the 
measurement by a few tenths of mK. We adopt the simulated 
values as a conservative estimate of the error due to misesti- 
mation of the detector offsets in Table [2] 

Although the polarized beam shape B has been measured to 
be symmetric between bolometers sharing a feed, the error of 
this measurement is approximately 2.5% of the beam width. 
The uncertainty in the individual detector beam widths could 
cause a systematic error in this analysis. As there is no ev- 
idence for a systematic shift in the beam widths, to quantify 
this error we have performed simulations where the pair beam 
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TABLE 2 

Cen A measurement uncertainty budget. 



Error on parameter 


EGlOOGHz 


S^lOOGHz 


S2l50GHz 


E^150GHz 


#100GHz 


#150GHz 


Beam pair offsets 


9mJy 


13 mJy 


lOmJy 


18mJy 


0°15 


0°22 


Differential beam shape 


4 mJy 


5 mJy 


6 mJy 


6 mJy 


0°06 


0°11 


Relative pair polarization angle 


9mJy 


6 mJy 


5 mJy 


4 mJy 


0°12 


0°09 


Background structure 


1 mJy 


OmJy 


4 mJy 


32mJy 


0°01 


0°29 


Source variability 


14 mJy 


31 mJy 


56 mJy 


68mJy 


0°28 


1°05 


Calibration angle 


35mJy 


35mJy 


28mJy 


28mJy 


0°50 


0°.50 


Total systematic 


39mJy 


49mJy 


64 mJy 


83mJy 


0°61 


1°23 


Random 


65 mJy 


57 mJy 


64 mJy 


70 mJy 


0°91 


1°17 


Total error 


76 mJy 


75 mJy 


91 mJy 


109 mJy 


1°09 


1°69 



widths are scattered about their known values using a Gaus- 
sian randomization with la = 2.5%. In addition, the relative 
detector pair polarization angle uncertainty is also a poten- 
tial source of polarization error. Simulations are performed 
assuming a la = 1° Gaussian scatter on the measured detec- 
tor pair relative angles, which is the estimated per pair uncer- 
tainty on ip (see H09). The results of both of these simulation 
sets are listed in Table 12 

As discussed in P09, we estimate the polarization leakage 
of our detectors to be e = 0.08 ± 0.015. Random errors in e 
will average down, while a systematic change merely shifts 
the absolute calibration by a factor of ~ 2a e . As the uncer- 
tainty on e is < 0.02 this is sub-dominant to the uncertainty on 
the overall absolute calibration. Sidelobe pickup of polarized 
emiss ion from the ground may also be a concern dBrown et alJ 
2009). However, both use of the fifth-order polynomial fil- 
ter and the small angular extent of the source suppresses the 
magnitude of this effect to the < 10 /iK level, so it is not a 
significant source of systematic error in this measurement. 

Table [2] summarizes the error budget for this measurement 
including all important systematic errors. The total systematic 
and total errors have been calculated using the square root of 
the quadrature sum of each of the individual errors. 

7. CONCLUSION 

We have presented measurements of the 4° x 2° region 
centered on the radio source Centaurus A with QUaD, a 
millimeter- wave polarimeter whose absolute polarization an- 
gle is known to ±0°.5. Systematic errors from astronomical 
sources such as the Cen A field's background structure and 
temporal variability from the source's nuclear region have 
been estimated. Simulations have been performed to assess 



the effect of misestimation of the instrumental parameters 
on the final measurement. After correcting for known in- 
strumental effects, we find that the total (g, U) of the inner 
lobe region is (1.00 ± 0.08,-1.72 ± 0.08) Jy at 100 GHz and 
(0.80±0.09,-1.40±0.11)Jy at 150 GHz, leading to polar- 
ization angles of -30?0 ± 1?1 and -29?1 ± 1?7. 

Future millimeter polarimeter s, including CMB experi- 
ments, will require a well characterized astronomical source 
against which to compare their instruments' laboratory deter- 
mined polarization properties. Cen A represents one of the 
best candidates for a stable, compact, highly polarized source, 
and the measurements presented here will allow calibration of 
such instruments to an uncertainty of ~ 1°. 
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